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Abstract 

We study the effects of chirality at the segment scale on the thermodynamics of block copolymer 
melts using self consistent field theory. In linear diblock melts where segments of one block prefer 
a twisted, or cholesteric, texture, we show that melt assembly is critically sensitive to the ratio of 
random coil size to the preferred pitch of cholesteric twist. For weakly-chiral melts (large pitch), 
mesophases remain achiral, while below a critical value of pitch, two mesocopically chiral phases 
are stable: an undulated lamellar phase; and a phase of hexagonally-ordered helices. We show that 
the non-linear sensitivity of meso-scale chiral order to preferred pitch derives specifically from the 
geometric and thermodynamic coupling of the helical mesodomain shape to the twisted packing of 
chiral segments within the core, giving rise to a second-order cylinder-to-helix transition. 
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Widely admired for their value to nanotechnology block copolymer (BCP) melts 
are also well-studied prototypes of soft-molecular assembly. Even melts of the simplest, 
diblock architecture exhibit a rich spectrum of periodic, long-range ordered states, as a con- 
sequence of the competition between flexible chain entropy and demixing of unlike chemical 
domains |2j. The nearly limitless flexibility of molecular design of copolymers, as well as a 
range of powerful experimental and theoretical methods for studying phase behavior, have 



'actors such as molecu- 
^ dictate properties of 



provided an extensive platform for testing how changes in many 
lar shape {3, 4] or the nature and range of inter- molecular forces 
long-range ordering in self-assembled systems. One aspect of self-assembly that is crucial 
in a broader class of systems from liquid crystalline materials to biological assem- 

blies js, 9] is the influence of molecular-scale chirality on long-range ordering. Despite the 
tremendous progress in understanding block copolymer thermodynamics in recent decades, 
the generic influence of chirality at the monomer scale on the meso-scale phase behavior is 
largely unknown. 

Indeed, recent experimental studies by Ho and coworkers provide direct evidence that 
chirality at the scale of monomers has a profound impact on the meso-scale assembly of 



certain block copolymers 



10Nl2|. Unlike the standard cylinder (C) phase of achiral copoly- 



mers, melts of polystyrene- 6-poly(L or D)lactide (PS-PLLA or PS-PDLA) assemble into 
hexagonally-ordered arrays of helices, denoted as the H* phase, where the handedness of 
helices is shown to switch with a reversal of monomer chirality 12J|. Even in this homochi- 
ral system, transfer of chirality from the molecular- to the meso-scale is not automatic: at 
lower compositions of the chiral block, the achiral C phase is also observed 11] . These 
observations beg a number of unanswered questions about the influence of chirality on the 
block copolymer phase diagram. Specifically, how does the mesoscopic packing of chains in 
the assembly couple to a thermodynamically preferred chiral packing in a way that reflects 
molecular handedness? What key parameters dictate the stability of meso-chiral assemblies? 
Finally, along with H*, what other meso-chiral phases are stable in phase diagram of chiral 
block copolymer melts? 

In this paper we investigate the equilibrium phase behavior of melt of diblock copolymers 
possessing a single chiral block. Our approach is based on a recent self-consistent field theory 
(SCFT) approach to block copolymer melts that considers the mean-field cost for gradients 
in vector order parameter, t(x), the local mean segment orientation in ordered phases [131 ] . 



2 



The chiral nature of the polymers is modeled as a thermodynamic preference for a cholesteric 
twist of the orientational order associated with the chiral block. We map the phase diagram 
of weakly-segregated melts in terms chiral segment composition /, degree of segregation 
and strength of chirality, as measured by go — 27r/p, where p is the preferred pitch of 
the cholesteric helix of chiral segments. Above a critical value of go, chirality dramatically 
alters the phase diagram leading to a stable H* phase, as well as new morphology, a phase 
of undulated layers, driven by cholesteric segment twist. Finally, we analyze the C-to-H* 
transition based on both our SCFT results and a strong-segregation picture of chain packing 
to show that non-linear sensitivity of chirality transfer to go derives specifically from the 
geometric interplay between mesodomain formation and inter-segment twist, giving rising 
to a second- order phase transition between these two phases. 

We consider incompressible melts of AB diblock copolymers possessing N Kuhn segments, 
f'N of which belong to the chiral A block. The total free energy of the melt (in units of 
ksT) is written as, 

F = p 0X J dVM*)Mx) ~ Schain + F\ (1) 

where incompressibility requires a (x) + 0b (x) = 1. The first term represents the mixing 
enthalpy of unlike segments, where 0a ( x ) and 0b ( x ) represent the local volume fractions of 
A and B monomers, \ is the Flory-Huggins parameter describing unlike segment interactions 
and Pq 1 is the volume per segment. The second contribution denotes the entropic free energy 
cost for polymeric chains to deviate from random-walk (Gaussian chain) statistics. While 
these terms are familiar ingredients to the standard Gaussian-chain model of copolymer 

nn 

melts [14|,H5|], the third term represents the chiral coupling between gradients in the segment 
orientation and chain conformation. Here, t(x) is a vector order parameter corresponding 
to the local average of t a , the orientation of the ath A-block segment in the melt, 

t(x) = p 1 ^t Q 5(x-x Q ), (2) 

a=A n 

where the sum is carried out over all A segments on all chains [13j. To describe the 
orientation-dependent interactions between the A segments in a coarse-grained way, we 
adopt a g eneralization of Frank elastic energy, standard to the theory of liquid crystalline 
order 16], describing square gradient costs for order-parameter variations in a cholesteric 

F* = ^ J ^{^(V • t) 2 + K 2 (V x t) 2 + 2q Q K 2 t ■ (V x t) + K' 2 [t • (V x t)] 2 }, (3) 
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FIG. 1: In (a), the phase diagram for fixed segregation strength, \N = 14. In (b), the phase 
diagram a fixed degree of chirality, q = 3.6. Phases are labeled as in the text (D = disordered 
melt). 

where Ki, K% and K' 2 are elastic constants. A lack of chiral symmetry at the segment scale 
corresponds to the case go 7^ 0, and a thermodynamic preference for non-zero cholesteric 
twist of segments, t ■ (V x t) ^ [6|. In particular, for K' 2 = the minimal free energy 
texture has pitch p = 27i/q . While the first line of eq. ([3]) includes the square-gradient 
terms allowed by symmetry to second order in t(x), it is necessary to include the higher 
order term [K' 2 ^ 0) for stability. 

The mean-field, SCFT analysis of copolymer melts possessing only density-dependent 
interactions is well-known [14], and for brevity, we only summarize the key additional ele- 
ments required by the self-consistent coupling to t(x) (see ref. [l3] for details). To determine 
composition and orientation profiles in the melt as well as the mean-field chain entropy in 
the ordered phases, segment-segment interactions are replaced by the self- consistent fields 
generated by their mean distributions. We define u>a(x) and »b(x) to denote the density- 
dependent fields acting on A and B blocks, respectively, and W(x) denotes the orientation- 
dependent field acting on the A block. A modified diffusion equation describes the random- 
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walk statistics on chain conformations subjected to the spatially varying fields, 

dq = i f [aV + W(x)]V^ A (x)g, *</ 
ds j fa 2 V 2 g- WB (x)g, s>/' 

where g(x, s) describes the statistical weight of a section of chain spanning from the free end 
of the A block (s = 0) to the (siV)th segment situated at x and a is the mean-segment size. 
A similar equation describes the weight, q'(n, s), of "diffusing" from the other free end of the 
chain (s = 1) to the same point, such that the probability of the (sA)th segment of an entire 
chain at x is g(x, s)q ! (yi, s)/Q, where Q = J dVq{^x, s)gt(x, s) is the single-chain partition 
function. The distributions g(x, s) and (^(x, s ) completely determine the mean composition 
profiles as well as the mean-orientational profile of the A-block segments according to 

Y a ff r -I 

*( x ) = 7^ ds s)Vg t (x, s) - g t (x,s)Vg(x,s) - 2tr 1 W(x)g(x, ^(x, s) . (5) 

6Q Jo i J 

The condition for mean-field in the melt requires that self-consistent fields are related to 
the mean composition and orientation profiles by ?«a(x) — u>b(x) = x[0b(x) — 0a(x)] and 
Wi(x) = po5 F* / 5ti(x) . These equations, along with the incompressibility condition, form 
a closed set whose solutions correspond to the equilibrium (mean- field) states of the melt. 
The free energy of mean- field solution can be calculated directly from 0a, <Pb arid t(x) and 
a relation for the entropy (per chain), Sc^n/n = \nQ + V' 1 f cIV[wa4>a + wb0b + W • t]. 

The self-consistent equations are solved numerically on a three-dimensional, periodic grid 
(of typical dimensions 20-35 on a side). The distribution functions, g(x, s) and ^(x, s), are 
solved by numerical integration of the "diffusion + drift" equations, allowing for a direct 
calculation of 0a, 0b and t profiles for given self-consistent field distributions. In turn, 
the self-consistency relations are solved for w\, wb and W via an iterative, over-damped 
relaxation algorithm [3]. To map the phase behavior, the minimal free-energy solution 
is identified for a given set of thermodynamic parameters. For this study, we consider 
melts described by the one Frank constant approximation, K\ — K 2 = NK' 2 = K, we set 
Na 2 K = 1/2 and we focus on the regime where the chiral block is the minority component, 
/ < 0.5. Three thermodynamic parameters remain to govern melt behavior: degree of 
segregation, xiV; chiral fraction of the chain, /; and the degree of chirality q = q Q N l / 2 a^ 
proportional to the ratio of (ideal) chain size to preferred cholesteric pitch of the chiral 
block. 
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FIG. 2: SFCT results for the UL* phase at \N = 14, / = 0.44 and q = 3.9: in (a), the 3D volume 
of chiral domain is shown; in (b), an edge-on view of the segment orientation within the chiral 
layers (arrows indicate t-field orientation). A cartoon schematic of the chain packing is shown in 
(c), with net polarization of chiral block segments shown as red arrows and achiral block shown as 
blue coil. 

In Fig. [T^, we show the phase behavior in the (/, q) plane for fixed, intermediate segre- 
gation strength %iV = 14, while Fig [lb presents the (/, xiV) plane for a fixed degree of 
chirality, q = 3.6. In the former plane, we find that the phase diagram is divided into two 
regimes, weakly-chiral melts, q < 3.45 and strongly-chiral melts, q > 3.45. For weakly-chiral 
melts, we predict the standard sequence of BCP phases with increasing composition: a BCC 
lattice of spherical domains (S), a hexagonal lattice of cylinders (C), a cubic double gyroid 
(G) and a lamellar phase (L). Despite preference for twist in the system, the phases of the 
weakly-chiral regime are not chiral on the mesoscale, and consequently, the phase boundaries 
separating them do not show dependence on q. In contrast, in the strongly-chiral melts at 
the highest composition range we find that the achiral C, G, L phases are overtaken by two 
phases that do not appear in achiral BCP phase diagram, helical cylinders H*, and an un- 
dulated lamellar phase UL*. From the (/, xN) section of phase space at q — 3.6 in Fig. [lb, 
we see that the stability of the mesochiral H* structure is not dictated by degree of chirality 
alone, as this morphology is found to have a window of thermodynamic stability only above 
a critical degree of segregation, xiV ~ 12.5. 

The thermodynamic stability of both the H* and UL* phases derives from the threading 
of cholesteric texture through the minority domains of either morphology, which is apparent 
upon inspection of the segment orientation distribution in these phases. In the UL* phase 
Fig. [21 t(x) twists around a cholesteric pitch axis running along the layer, this direction also 
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FIG. 3: SFCT results for the H* phase at %A = 14, / = 0.345 and q = 3.6: in (a), the 3D volume 
of chiral domain is shown; in (c) segment orientation within the consecutive cross-sectional slices 
through the chiral domain shown in (b). In (d), a schematic depiction of coupling of segment 
orientation in cross section and the bending of twist of the backbone of the tubular core domain. 
As in Fig. [21 the net polarization of the core block segments is depicted as a red arrow. 

being the direction of lamellar layer undulation. Unlike the standard "double-layer" lamellar 
morphology, the UL* phase is composed of single layers whose chiral blocks are extended 
and "highly-polarized" . The density of the achiral blocks emerging from either side of these 
layers oscillates due to the rotation of layer polarization, which in turn drives a periodic 
bending, or undulation, of lamella whose periodicity is locked to the cholesteric twist. In 
H* (Fig. [3]), we find the segments of minor (chiral) domain to be polarized along direction 
normal to the bending of the tubular domain, orienting towards the outside of the helix. As 
this direction rotates along the pitch axis of the helix, overlaying consecutive sections of the 
helical core show a cholesteric twist that threads along the backbone of the helical domain 
(Fig. Eh). 

SCFT results demonstrate that chirality at the segment scale is a necessary, but not 
sufficient condition, for equilibrium stability of mesoschiral phases like H*, consistent with 
experimental observations [11]. To explain the critical sensitivity of chirality transfer on 
/, \N and q we consider a simplified analysis of the stability of the H* phase over C, the 
achiral cylinder phase. At the heart of the chirality transfer from segments to mesoscale in 
the H* phase is the geometric coupling of bending of the tubular domain and polarization of 



segment orientation in the cross section, shown schematically in Fig. [3jd. Because the radial 
segment orientation within the core of the C mesodomain is isotropic — orienting equally in all 
directions in the planar section — a simple rotation, or twist, of consecutive sections generates 
no cholesteric twist, as measured by t ■ (V x t). Bending of the tubular domain breaks the 
isotropic symmetry of the C domain, tending to orient segments along the unit normal of 
the center line of the tube. We adopt a strong-segregation perspective [r3] to estimate the 
degree of orientation by dividing the core domain into a series of wedges containing A-block 
chains extending radially from the centerline to the AB interface along the vector Ra, and 
note that the mean segment orientation within such as wedge is Ti^/(fNa). Accounting 
for the excess volume of wedges on the outer side of the bent tubular domain, we average 
over all wedge orientations to find mean segment orientation t ~ —a(f)K>hR 2 /(Na), where 
k and n are the curvature and normal of the tube center line, R is the outer domain radius 
of the tube, and a(f) = — l)/6 ~ f~ 1 ^ 2 /2. Hence, the helical rotation the normal 

direction along the backbone of the tubular domain, leads to a cholesteric twisting of the 
polarized segments in the core, t • (V x t) ~ — |t| 2 r, where r is the torsion of helix and we 
can estimate the chiral free energy gain (per chain) of the H* phase by replacing t(x) with 
its cross-sectional mean, 

where we set r = go to maximize chiral free-energy gain. Counteracting the tendency to 
polarize then twist the cores is the mechanical cost of bending the cylindrical microdomains 
(per unit cylinder length), F^ cn d/L = (B/2)k 2 . We estimate the modulus from the de- 
formation cost of lamellar microdomains [18| as B = /3(/)p ~ 1 A(xA^) 1 / 3 i? 4 , where /3(f) is 
a numerical prefactor which we assume to be weakly dependent on /. Combining these 
estimates we find the free-energy dependence of the H* phase in the small curvature limit, 

F(H*)/n^c(q 2 c -q 2 )K 2 + O(K*), (7) 

where c ~ KR A and Kq 2 m (^iV) -1 / 3 assuming the strong-stretching domain size R ~ 
(xNy^N 1 / 2 !!. Hence, the mechanical cost and chiral gain of helical buckling of core domain 
both grow with square curvature. While for weakly-chiral melts, go < Qc, the straight 
(k = 0) and achrial domains of the C are stable, while for g > g c the core becomes 
thermodynamically unstable to helical bending. Near to the critical degree of chirality we 
expect curvature of the core domain to grow continuously from zero as k ~ |go — gd 1 ^ 2 - Due 
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FIG. 4: Volume-averaged cholesteric twist of SCFT results at %iV = 14, showing a phase transition 
from C to H* with q and increase of q c with decreasing chiral block composition. 

to the geometric coupling between segment orientation and domain bending, we find that 
chirality transfer is highly non-linear process in tubular phases of copolymer melts, taking 
place as a second-order phase transition from the C-to-H* phase. The signature of critical 
sensitivity to go is shown in the volume average of cholesteric twist extracted from the SCFT, 
shown in Fig HI Consistent with the thermodynamic argument above, we find no twist for 
go < g C) while above the critical degree of chirality SCFT shows linear increase in cholesteric 
twist, following the critical scaling t • (V x t) ~ k 2 ~ |g — g c |. 

In summary, SCFT studies of chiral block copolymer melts reveal the equilibrium phase 
behavior of melts is critically sensitive to the thermodynamic preference for cholesteric twist 
of the chiral domain. A drive for twist at the scale of chain segments indeed reshapes the 
structure and symmetry of the ordered phases on much larger length scales. Importantly, 
the experimentally observed H* morphology is determined to be stable in the equilibrium 
phase diagram of chiral melts, and its thermodynamic stability derives directly from the 
geometric coupling of domain distortions and inter-segment twist. While further studies are 
needed to map its full influence, there is little question that molecular chirality has a direct 
and potent impact on the structure and thermodynamics of block copolymer melts. 

The authors are grateful to H. Johnston for extensive and invaluable discussions. This 
work was supported by the MRSEC on Polymers at the UMass (WZ, GG) and DOE-DE- 
FG02-45612 (TPR). 
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